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"-^ . In preparation for a set of hydrodynamical simulations of core-collapse supernovae and protoneu- 

OO ' tron star winds, we investigate the rates of production and thermalization of i/^ and Vr neutrinos in 

Cn , dense nuclear matter. Included are contributions from electron scattering, electron-positron annihilation, 

nucleon-nucleon bremsstrahlung, and nucleon scattering. We find that nucleon scattering dominates 
electron scattering as a thermalization process at neutrino energies greater than ~15 MeV. In addition, 
O^ ' nucleon-nucleon bremsstrahlung dominates electron-positron annihilation as a production mechanism at 

^^ , low neutrino energies, near and below the f^ and i>r neutrinosphere. 

xj ' Furthermore, we have begun a study of steady-state general relativistic protoneutron star winds em- 

ploying simple neutrino heating and cooling terms. From this analysis we obtain acceleration profiles as 
well as asymptotic lepton fractions and baryon entropies essential in assessing the wind as a potential 
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o . 

^^ , site for r-process nucleosynthesis 

'^ [ 1. Introduction 

Q ' The cores of protoneutron stars and core-collapse supernovae are unique environments in nature. They 

M I are characterized by mass densities of order ~ lO'^^ — 10^^ g cm~'^ and temperatures that range from 

t/^ ■ ^ 1 to 50 MeV. At these temperatures and densities neutrinos of all species are produced in proliferation 

via electron-positron annihilation (e+e^ ^-> i/p), nucleon-nucleon bremsstrahlung, and plasmon decay 
ilpi *~* ^i^)- While these processes contribute for the electron types (z^eS and Ves), for them the charged- 
current absorption and emission processes VeU <-» pe~ and DeP ^^ ne^ dominate. Neutrino-electron 
; ^ ■ and neutrino-nucleon scattering also contribute to the total opacity. All of these interactions combine 

^ I to couple the neutrinos to dense nuclear matter, affecting energy transport from the core where the 

neutrinos are diffusive to the more tenuous outer layers where the neutrinos begin to free-stream. Indeed, 
the neutrino heating in the semi-transparent region behind the shock is now thought to be an essential 
ingredient in igniting the supernova explosion itself [pm[. Furthermore, a neutrino-driven protoneutron 
star wind is thought to be a general feature of the core-collapse phenomenon and has been proposed as 
a site for r-process nucleosynthesis. 

Essential in understanding the mechanism of Type-II supernovae is an accurate prediction of the 
production spectrum for each neutrino species. To this end, in §2 we present results of a thermalization 
and equilibration study of fj, and r neutrinos in dense matter. 

Separately, in §3 we present preliminary results from calculations of steady-state protoneutron star 
winds in general relativity. We include velocity profiles, asymptotic entropies, and expansion timescales 
central in assessing this site as a candidate for r-process nucleosynthesis. 

2. Thermalization and Production 

For f^, and I'r neutrino types (collectively V^s'), which carry away 50—60% of the ^^ 2 — 3 x 10^'^ ergs 
liberated during collapse and explosion, the prevailing opacity and production sources are i^^-electron 
scattering, i/^-nucleon scattering, e~^e~ annihilation, and nucleon-nucleon bremsstrahlung. The charged- 
current reactions dominate the electron-type transport so completely that we do not consider them here. 
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Our focus is on (1) the role of i/^-nucleon scattering relative to i^^-electron scattering in thermalizing z/^ 
neutrinos and (2) the importance of bremsstrahlung as compared with e^e~ annihilation as a producer 
of Vfii^fi pairs. 

Supernova theorists had long held ^ that i^^-nucleon scattering did not aid in thermalizing any neutrino 
species. While the process was included as a source of opacity |^g it served only to redistribute the 
neutrinos in space, not in energy. In contrast, i/^-electron scattering was thought to dominate as a 
thernializer of v^s. In addition, the only Uf^D^ pair production mechanisms employed in full supernova 
calculations were e^e^ <-> I'^f'^ and plasmon decay (7pi <-> v^v^) H. Nucleon-nucleon bremsstrahlung, 
while recognized as a late-time cooling mechanism for more mature neutron stars |qM| , was neglected in 
supernova theory. Recent developments, however, call both these practices into question and motivate 
a re-evaluation of these processes in the supernova context. In the last few years, analytic formulae for 
I'^-nuclcon have been derived that include the full kinematics and final-state Pauli blocking, at arbitrary 
nucleon degeneracy, at the temperatures and densities encountered in the core of a core-collapse supernova 
l8l-0| . These efforts reveal that the average energy transfer in i/^-nucleon scattering may surpass previous 
estimates by an order of magnitude and, hence, that this process may compete with j/^-electron scattering 
as an equilibration mechanism |T2-n3] . Similarly, estimates of the nucleon-nucleon bremsstrahlung rate 
have been obtained P,0,Qq,O| which indicate that this process can compete with e+e^ annihilation in 
the dense core. 

In order to compare these scattering and production processes directly, we solve the Boltzmann equation 
for the time evolution of the neutrino phase-space distribution function {J-^) in an idealized system with 
no spatial or angular gradients. We consider an isotropic homogeneous thermal bath of scatterers and 
absorbers held at constant temperature, density, and electron fraction. For the scattering processes, 
we begin the equilibration calculation at t = with a v^ distribution function with a characteristic 
temperature of twice that of the surrounding matter. We then evolve this distribution function using 
the full collision integral of the Boltzmann equation, with the structure function formalism of Burrows 
and Sawyer (1998) [|| and Reddy et al. (1998) g. When equilibrium is reached the final distribution 
is Fermi-Dirac, at the temperature of the surrounding matter, with chemical potential set by the initial 
neutrino number density (n^), conserved to better than 0.0001% throughout the calculation. For the 
production and absorption processes we begin with zero phase-space occupancy for both j/^ and i/^ 
at all energies and let bremsstrahlung and e^e~ annihilation each build an equilibrium distribution of 
neutrinos and anti-neutrinos. e^e~ annihilation is calculated in a Legendre polynomial expansion [g|. 
The production rate via nucleon-nucleon bremsstrahlung is calculated in a one-pion exchange model with 
arbitrary nucleon degeneracy ||l8|,|9[. As a check to the calculation, the final distribution should have a 
characteristic temperature of the ambient matter with zero neutrino chemical potential. 

The left panel of Fig. n^ shows the thermalization rates for i/^n and v^e~ scattering, defined in terms 
of the average energy transfer {to) at that energy, in equilibrium. The calculation was performed with 
T ~ 6.1 MeV and p ~ 1.1 x 10^^ g cm~^. These thermodynamic conditions are representative of the v^ 
neutrinosphere, the semi-transparent regime where the neutrinos begin to decouple from the matter and 
free-stream to infinity. Most noticeable in this graph is the fact that while v^e~ scattering dominates 
at low energies (< 10 MeV), at modest and high energies v^n scattering competes with or dominates 
thermalization. Throughout our calculations, at a variety of densities, temperatures, and compositions, 
we find this behavior to be generic for the two scattering processes. Indeed, the point where v^n scattering 
begins to overwhelm v^e~ scattering seems always to fall between approximately 10 and 20 MeV. 

The right panel of Fig. Q reveals the same type of systematics for the production and absorption pro- 
cesses. In this case, however, we plot the total differential emissivity for nucleon-nucleon bremsstrahlung 
and e^e~ annihilation in equilibrium for two different thermodynamic points taken from a one-dimensional 
core-collapse simulation |3|Jl9(| . Bremsstrahlung is clearly the dominant production mechanism at low neu- 
trino energies. In the interior, at densities of order 10"'^'^ g cm~'^, e+e~ annihilation begins to compete 
only at neutrino energies above 50 MeV. In our time-dependent calculations, we find that bremsstrahlung 
always dominates production below 10-20 MeV at all points in a representative collapse profile. 

3. Neutrino-Driven Protoneutron Star Winds 

A complete and self-consistent theory of the origin of all the elements has been the grand program 
of nuclear astrophysics since the field's inception. The r-process, or rapid neutron capture process, is 



a mechanism for nucleosynthesis by which seed nuclei neutron capture on timescales shorter than those 
for (3~ decay. With a sufficient neutron flux, capture continues to very neutron-rich isotopes and to 
the heaviest elements (e.g., Eu, Dy, Th, and U) producing unique abundance peaks aX A ^ 80, f30, 
and 195 po|-p^. The r-process is only quelled when photodisintegration timescales approach those for 
neutron capture. After the intense neutron flux lessens, (3" decay populates the primary stable isobar for 
a given atomic number. While the relevant nuclear physics is fairly well understood, the astrophysical 
site, which must exist in order to produce the elemental abundances we find in nature, is not known. The 
viability of a site for r-process nucleosynthesis hinges on three characteristics: the asymptotic entropy 
per baryon (s/), the electron fraction (Kg), and the dynamical timescale (rdyn). One proposed site is the 
protoncutron star wind that emerges after core collapse and shock reheating during a supernova P,E3,e3 . 
Both numerical and analytic studies of the conditions in this neutrino-driven wind have been carried 
out previously [p5 26 . Early calculations based on realistic supernova models produced interesting nu- 



cleosynthesis and appreciable r-process yields, but overproduced nuclei near A^ = 50 (^^Sr, ^^Y, and 
^"Zr) |2jJ2q]. This problem was overcome by fine-tuning Ye in these simulations p3], but no consensus 
on the other parameters (particularly, Sf) has yet been reached. The pioneering analytic work of ref. 
p6[ using simple wind models showed that entropies fell short by a factor of ~ 2 — 3 of those needed 
for the dynamic timescales and lepton fractions achieved. In other studies s/ was artificially enhanced 
by a factor of ^^ 5 to achieve proper solar r-process abundances [ESi. The fact that these calculations 
indicate an Sf too low for r-process nucleosynthesis does not exclude the wind as a potential site. A slight 
increase in energy deposition after the wind's acceleration phase p6[ or the effects of general relativity 
have been shown to decrease Tdyn and increase s/ [Q, both favorable to nucleosynthesis. Recent general 
relativistic steady-state and hydrodynamical studies indicate that winds can generate all three r-process 
abundance peaks only when the protoncutron star is quite massive ~ 2 Mq |31| , |32| and the total neutrino 
luminosity large (~ 10^^ erg s~^). These condition produce modest entropies (~ 130 baryon"^ kg^), but 
short Tdyn ^ 6 ms. 

The wind equations can be reduced to three ordinary coupled, critical differential equations for the 
evolution of temperature (T), mass density (p), and velocity [v) in radius. Yg is held constant. Solving 
these equations constitutes an eigenvalue problem. The eigenvalue sought is the mass outflow rate 
M ■= Airr'^pv, or, alternatively, the critical radius {Re) where v{Rc) = Cs, the local speed of sound. 
In practice, we impose two boundary conditions at the protoncutron star surface (Ro), which we take to 
be the Ve neutrinosphere: (1) T(Ro) = T^^ and (2) Tv{Ro) = — J K,^pdr = 2/3, where T,y is the neutrino 
optical depth. A third constraint must be imposed if we are to close the system of equations. This we take 
to be the critical condition, v{Rc) — Cs, which defines the outer radial boundary. The system, now well- 
posed, can be solved using a relaxation technique on an adaptive radial grid [ p3[ . The code then adjusts 
the radial mesh in a Newton-Raphson sense in order to fulfill all boundary conditions simultaneously. 
Once the critical point is determined, we use I'Hospital's rule to bridge it and then integrate to infinity 
using a simple Runga-Kutta scheme. 

The results of a preliminary and representative general relativistic calculation are shown in Fig. 0. The 
left panel shows velocity profiles for a variety of Pg luminosities, for a protoneutron star mass M = IAMq 
and radius Ro = 10 km. This range of neutrino luminosities is indicative of the first '^ 8 — 10 seconds of 
the protoneutron star's life. The right panel is a plot of the asymptotic entropy per baryon per Boltzmann 
constant versus Tdyn for the same protoneutron star. Note that these calculations were carried out with 
constant Y^ = 0.302. 

In general, we find that higher entropies and shorter dynamical timescales result from the use of general 
relativity instead of Newtonian gravity. For the 1.4 Mq object we consider here, however, we do not 
achieve entropies and timescales appropriate for r-process nucleosynthesis |29, 3l] , |3^ . 



4. Summary and Conclusions 

Our results for equilibration via f^-electron scattering and i/^-nucleon scattering demonstrate that 
the latter competes with or dominates the former as a thermalizer for neutrino energies > 10 MeV 
for p > I X 10^^ g cm^'^ at all temperatures. At neutrino energies > 30 MeV, the difference at all 
densities and temperatures is approximately an order of magnitude. For the production and absorption 
processes, we find that nucleon-nucleon bremsstrahlung, at the average energy of an equilibrium Fermi- 
Dirac distribution at the local temperature, is 2 orders of magnitude faster than e"'"e~ annihilation at 



T '^ 15 MeV and p ~ 10"^^ g cm"'^. Only for p ~ 10^^ g cm"'^ and T ~ 6 MeV does e+e" ^-> v^D^ begin 
to compete with bremsstrahlung at all energies. We conclude from this study that the emergent i'^ and 
Vt spectrum is (1) brighter and (2) softer than previously estimated |16[| . The former results from the 
inclusion of the new pair emission process, nucleon-nucleon bremsstrahlung. The latter is a consequence of 
both the increased energy coupling between the nuclear and neutrino fluids through i/^-nucleon scattering 
and the fact that bremsstrahlung dominates e+e~ annihilation near the neutrinospheres at the lowest 
neutrino energies. 

In addition, our first step in the calculation of realistic protoneutron star wind models has been success- 
ful; wc have created a robust technique for solving the steady-state eigenvalue problem and have confirmed 
the results of other researchers qualitatively. We plan to explore the parameter space of protoneutron 
star winds exhaustively with an eye toward implementing full neutrino transport in a hydrodynamic 
simulation of wind emergence and evolution. 
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Figure 1. Left panel: the thermalization rate T (s~^) for both j/^-electron and I'^i-neutron scattering as 
a function of neutrino energy Ei, (MeV) in equilibrium. Right panel: the differential emissivity versus 
El, for both nucleon-nucleon bremsstrahlung and e+e^ annihilation for two thermodynamic points taken 
from a representative one-dimension core-collapse simulation 0] . 
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Figure 2. Left panel: general relativistie protoneutron star wind velocity v (cm s ^) profiles for six v^. 
luminosities: (A) 4.5, (B) 3.5, (C) 2.5, (D) 1.5, and (E) 1 x 10^^ erg s~^ The mass and radius of the 
central neutron star are M = 1.4Mq and Ro = 10 km, respectively. The electron fraction (Yg) is held 
constant at 0.302. The dots denote the critical point where the velocity of the wind is equal to the local 
speed of sound. Right panel: asymptotic entropy (s/) versus dynamical timescale (rdyn) for constant 
protoneutron star mass and radius, but as a function of v^ luminosity. 



